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Abstract 



A new model of self-organized criticality is proposed. An algebra of operators 
is introduced which is similar to that used for the Abelian sandpile model. 
The structure of the configurational space is determined and the number of 
recurrent states is found. 



Since the introduction of self-organized criticality (SOC) JT|, there has 
been considerable interest in the study of cellular automata which demon- 
strate how power-law correlations emerge during the evolution of extended 
dissipative systems. Of them, the Abelian sandpile model expresses most 
clearly the idea of dissipative dynamics, where a small disturbance exceeding 
a threshold grows and propagates through the system as an avalanche. 

Avalanches seem to be crucial for SOC, while the significance of a thresh- 
old is not completely clear. Indeed, some of the characteristics of sandpiles 
are purely diffusive even though one might expect more complicated be- 
haviour due to the presence of thresholds of stability. Thresholds in the 
forest fire model || are hidden in two parameters which are the probabilities 
of growth and ignitions. In the Bak-Sneppen model M, a threshold value 
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appears as a result of infinitely long evolution. 

In this Letter, I propose a cellular automaton model exhibiting SOC and 
containing no threshold parameters. However, the model has a common alge- 
braic structure with the Abelian sandpile and, therefore, a similar structure 
of recurrent configurations corresponding to the critical state. The parallel 
consideration of these models should elucidate the critical dynamics of both. 

Consider a two dimensional square lattice C of size L x L. Each site % of 
C is characterized by the radius-vector r« with integer Cartesian coordinates 
(xi,yi) and by one of the unit vectors ej(l), ej(2), e,(3), ej(4) directed up, 
right, down or left from % . At each discrete moment of time, one drops a 
particle on a site of C chosen at random and allows it to walk by the following 
rules: 

(i) at each step, the particle coming to a site j ; £ C turns the 
vector Gj{y) clockwise by the right angle: 

where \x = v + l(modA). 

(ii) performs the unit step in the direction e,-(/i) to the neigh- 
bour site j : 

r f = r j + e j(v) 

(iii) if the walk reaches a boundary site j and the new position 
j' is outside the lattice, the particle leaves the system. 

The walk is assumed to be quick enough to be completed by the next 
discrete moment of time. As a result of the walk, a configuration of vectors 
C specified by a unit vector e^u) on each lattice site i £ C transforms into 
a new configuration C which can generally differ from C by directions of 
vectors on sites visited by the particle. 

To describe the transformation resulting from dropping a particle on the 
site i , we define the operator acting on C and producing C : 

a t C = d (1) 

Theorem 1 The operator Oj exists, that is, for any C the walk started 
at i never enters a non-trivial cycle. 
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Proof: If a cycle contains a boundary site, the walk visits this site more 
than four times. Since one of these visits has as a consequence a step outside 
the lattice, no cycle can contain a boundary site. Next, consider a site one 
step away from the boundary. It cannot be visited in a cycle as, in this case, a 
particle must hit one of the the boundary sites. By induction, there can be no 
cycle containing an arbitrary site of the lattice and, therefore, the operator 
dj exists. 

The operators aj all commute. This property enables one to construct 
an Abelian group quite similar to that defined by Dhar for the sandpile 
cellular automata. 

Theorem 2 For arbitrary sites i and j and for any configuration of 
vectors C 

didjC = djdiC (2) 

Proof: Consider the updating procedure a^C = C' as a sequence of ele- 
mentary steps aj n ...aj 2 aj 1 with jx = i. The operator aij corresponds to the 
rotation of the vector ej(u) — *> Bj(n),n = u+l(mod4) and a consequent single 
step in the direction ej(/z). Then ciidjC can be written in the form 

where subscripts (1) and (2) refer to different particles. If j ^ j' , the opera- 
tors af and ot\, commute. If j = j' , they also commute due to identity of 
particles. Therefore, Oj and a,j commute. 

The commutativity rule enables us to consider several walks simultane- 
ously updating them concurrently. 

As usual, in the theory of Markov chains, we divide the set of all con- 
figurations {C} into two subsets, recurrent and transient. The first subset 
denoted by {R} includes those configurations which can be obtained from 
an arbitrary configuration by a sequential action by operators a. t . It follows 
from the definition that the subset {R} is closed under multiple action by 
operators Oj. Once the system gets into {R}, it never gets out under sub- 
sequent evolution. All nonrecurrent configurations are called transient and 
form the subset {T} which is the complement of the set {R}. 

Similarly to the avalanche operators of the sandpile model, the operator 
Oj has a unique inverse. The proof of this statement can be carried out in a 
way very close to the sandpile construction 0. 
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Theorem 3 For any recurrent configuration C G {R}, there exists a 
unique 

{afC) G {R} 

such that 

ai {ai l )C = C (3) 

Proof: Consider a standard recurrent configuration C* which can be chosen 
as the set of parallel vectors ej(l) on all % G C. 

Construct an identity operator EC* = C*. To this end, we take the 
product of boundary operators (i G B) multiplied by the product of corner 
operators <ij (i G A) 

n n °* (4) 

ieB ieA 

where B is the set of all boundary sites and A is the set of corner sites. 
Evidently, the operator (f|) does not change the standard configuration since 
it is nothing but a sequence of four successive rotations of arrows in all rows 
and columns of the lattice. Besides (f|), we can find another representation 
of E , having a nonzero degree of at any i G L. 

For instance, let us construct an identity operator E having operators 
in the next neighbours of the boundary sites. Replace a, by a,i and consider 
the product 

(nn«*y (5) 

\ieBieA / 

Let rij,z G C be the occupation numbers of particles resulting from C* after 
the action by the operator (J^).Each site a step away from the boundary 
receives at least one particle. Thus, the operator 

IK ( 6 ) 

i 

is the identity operator for C* having rii > at each next neighbour of 
edges. Repeating this procedure, we can construct the identity operator 
E(jii) having rij > at an arbitrary chosen site i £ £. 

Now, drop a particle on the configuration C* to obtain a new recursive 
configuration ajC*. The operator 

Pi = E( ni - 1) (7) 



4 



being combined with ai gives the identity operator again: 

p lCLl C* = C* (8) 

By definition, any recurrent configuration C can be obtained from C*\ 

C = P 2 C* (9) 

where P 2 is a product of the cij. Using the commutativity property, we have 

C = P 2 Pia i C* = aiP 2 P x C* (10) 

The intermediate result P 2 P\C* is the seeking configuration (a~ C). 

The proof of uniqueness doesn't differ from that for sandpiles [H]. Re- 



peating the construction (|T0l) , we can find a sequence C n £ {/?} such that 

(a,) B C B = C (11) 

Since the total number of configurations doesn't exceed 4 i2 , this sequence 
must enter a loop of length m > 1. This loop must contain C as C £ {/?} 
is attainable from an arbitrary point of the loop. We have a™C = C. Then 
a™ _1 C = a~ 1 C is the unique inverse. 

As all recurrent configurations can be obtained from an arbitrary one by 
the successive acting by operators c^, we can represent any C £ {R} in the 
form 

c = Y[( ai rc* (12) 

i<=C 

The L 2 -dimensional vector n labels the recurrent configurations. We can 
note that the operator a\ returns the vector &i{v) to the former position and 
initiates a motion of four particles at neighbouring sites of i. Therefore, the 
operator afaJ^aJ 2 aJ*aJ A , where ji,j 2 ,j3,j4 are the nearest neighbours of i, 
doesn't change the initial configuration and is actually the identity operator 
E. Using the Laplacian matrix A with elements 

Aij = 4 % = j 

A it j = — 1 i,j nearest neighbours 
Ajj = otherwise 
one can write down the identity operator in the form 

E = Y[af^ (13) 
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Eq. ( jl3D shows that two vectors n and n' label the same configuration of 
arrows if the difference between them is J2j m j^ij where rrij are integers. The 
L 2 - dimensional space {n} has a periodic structure with the elementary cell 
of the form of a hyper- parallelepiped with the base edges A^-, j = 1, 2, L 2 . 
Thus, the number of non-equivalent recurrent configurations is 

N = detA (14) 

which is Kirhhoff's formula for spanning trees || and Dhar's formula for 
sandpiles 0. 

The similarity to trees is not accidental. By construction, vectors involved 
into rotation don't form closed cycles if the motion is over. Due to finiteness 
of the lattice, each vector takes part in the motion some time or other. Each 
time the updating procedure is over, the collection of vectors reproduces the 
set of bonds of a tree. Thus, starting from an arbitrary noncorrelated set 
of vectors, we come to the set of spanning trees that are characterized by 
power-law correlations between different sites. 

The correspondence with sandpiles is not surprising as well. The algebra 
of the operators a, completely coincides with that of avalanche operators of 
the Abelian sandpile model 0. Moreover, the identity operator (|T3|) has the 
same form for both the models . This is the reason why the numbers of 
recurrent configurations coincide. 

Continuing the analogy between self-directing walks and sandpiles one 
can find the expected number Gij of full rotations of the vector at site j, due 
to the particle dropped at i 0. During the walk, the expected number of 
steps outside j is AjjGij whereas — J2kjtj GikAkj is the average flux into j. 
Equating both fluxes one gets 

^2G ik A kj = Sij (15) 
k 

or 

G« = [A- 1 ],, (16) 

The close analogy with sandpiles calls for a definition of avalanches in our 
model. The first step after landing of the particle generally leads to emegence 
of a cyclic configuration of arrows. As a result, the system leaves the recurrent 
set . It is natural to define the avalanche as a process of restoration of the 
recurrent state. It corresponds to successive rotations from the beginning of 
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the motion up to the moment when an acyclic configuration is restored for 
the first time and the structure of the spanning tree is reconstructed. The 
number of steps n is the duration of an avalanche, the number of different 
visited sites s is its size. If the first step doesn't lead to a cyclic configuration, 
we put n = 1, s — 1. 

When a closed loop appears on the given lattice, a branch of the dual 
tree gets disconnected on the dual lattice. The probability distribution P(s) 
of disconnected clusters follows the power law |7j 



where s is the number of lattice sites belonging to a cluster. Since the number 
of steps n which are neccessary to restore the tree is proportional to s, one 
can expect the similar power law for the avalanche distribution P{n). 

The proposed model admits a natural generalization on an arbitrary 
graph and arbitrary order of numeration of unit vectors directed to nearest 
neighbours of a given site. The main result flli]) remains unchanged where 
A should be defined as the Laplacian matrix of the given graph. 
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